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I. INTRODUCTION 



Although the field of quantum chaos found its first experimental challenges in atomic physics, during the past 
eight years several systems in condensed matter physics have been identified and studied from this point of view. In 
all of these cases the low-temperature electronic conduction properties of semiconducting heterostructures have been 
1 measured and analyzed to infer properties of the electronic spectrum and wavefunctions which may reveal quantum 
t-H , manifestations of chaos, or of the transition to chaos. Of these experiments, those which have involved patterning of 
two-dimensional electron gas systems to create more complex geometric scattering have in general not been amenable 
to detailed semiclassical analysis because the precise nature of the underlying confinement potential is not known. In 
contrast, the series of experiments performed at Nottingham Universityu and Bell LaboratoriesB which we will review 
and analyze here, involve tunneling through a quantum well in high magnetic field, for which the classical dynamics 
within the well can be accurately modelled and related to experimental observables. Moreover this system is unique 
in that three in situ experimental parameters, magnetic field, B, bias voltage V, and the angle between the electric 
and magnetic field, can be continuously varied so as to generate a transition from integrable to (essentially) fully 
chaotic dynamics. As such the system provides an ideal testing ground for semiclassical methods; we will find that 
is also raises some new questions in non-linear dynamics of some general theoretical interest. We will review and 
summarize below our theoretical study of this systemETEl, while also extending our earlier semiclassical treatment of 
I | the problemE 

We begin with some essential background. It is possible to fabricate atomically-flat layered materials, typically 
of the semiconductor GaAs and the insulator AlGaAs, such that the band-gap varies in a controlled manner along 
the longitudinal (growth) direction. This creates an effective one-dimensional potential for the electrons along this 
direction which confines the free electrons or holes (which are introduced by doping) in the semiconducting layers, 
known as "quantum wells". In the experiments in question, a thin "emitter" layer is fabricated with only one 
longitudinal quantum level filled by doping, followed by a barrier, a wide (~ 120 nm ) empty quantum well, a second 
barrier, and a conducting "collector" region. Thus one has in the longitudinal direction the potential profile of a 
"resonant tunneling diode" . As the voltage across the wide well is varied, states in the well can become aligned with 
the emitter level, and resonances in the tunneling current through the double-barrier structure will be measured if 
emitter and well states are coupled. In the experiments we will discuss, a large magnetic field is imposed, initially 
perpendicular to the layers, and has the effect of quantizing the electronic states within the layers into highly degenerate 
Landau levels. The electronic density in the emitter is such that typically only the lowest (n=0) Landau level is full, 
and so an exact selection rule forbids tunneling through any but the n=0 states of the well. Since both the emitter 
and well states have the same linear energy variation with B, one measures a series of peaks in the I-V characteristic 
of the device with voltage spacing corresponding to the longitudinal quantum states of the well, independent of the 
value of the magnetic field. Resonant tunneling diodes (RTDs) of this type have been well-studied for more than a 
. decade. 

The novel features of the more recent experiments arose when the magnetic field was tilted with respect to the 
electric field (tunneling) direction. We choose coordinates such that the normal direction to the layers is the z- 
O ■ direction, and the tilt angle is in the y-z plane, so that B = cos^z + smOy (see Fig. |l|). It was observed by Fromhold 
et alia that as the tilt angle was varied at high magnetic field (B=10T), the frequency of resonance peaks in the 
I-V would abruptly change at certain values of the voltage (typically this frequency approximately doubled). These 
workers initially attributed this frequency-doubling to a density of states (DOS) effect within the well. They showed 
numerically that at large tilt angles the electron dynamics in the well is almost completely chaotic, and that there exist 
new unstable periodic orbits which collide with both the emitter and collector barriers. These orbits were supposed 
to give rise (at certain voltages) to additional DOS oscillations via Gutzwiller's Trace Formula. A later numerical 
analysis of the quantum states within the well found that, unlike the usual level-clustering expected from the Trace 
Formula, in this case individual levels were scarred by the relevant periodic orbits once within each period of the 
Gutzwiller oscillation. Moreover, upon solving for the tunneling current by a quantum transfer matrix-method it was 
found that in many cases these scarred levels carried essentially all of the resonant tunneling currentou. 

Additional interesting features of this system were found when Muller et ain at Bell Labs undertook a careful 
study of the entire parameter space of magnetic field (0 < B < 12T), voltage ( < U < IV) and tilt angle, focusing 
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particularly on the transition region to chaos. They analyzed this large data set by plotting the peak locations in 
the B-V plane for fixed tilt angles between eleven and forty-five degrees. This representation of the data revealed a 
regions of peak-doubling and tripling which evolved in a complex manner with increasing tilt angle. These experiments 
motivated Shepelyansky and StonoS to analyze the global phase-space structure of electron motion within the well. 
Using a simplified model which neglected collisions with the emitter barrier they were able to estimate via the 
Chirikov resonance overlap criterion!!! the transition to global chaos in this system as a function of B,V,#. They also 
suggested that the onset of peak-doubling in the system was due to bifurcations of the simplest "traversing orbit" , 
an interpretation adopted and amplified upon by Muller et al.EI. However the classical analysis of Shepelyansky and 
Stone was not suitable for a developing a quantitative semi-classical theory, since it did not describe orbits which 
collide with both harpers, which dominate the tunneling spectrum according to both semiclassical theorya and exact 
quantum analysesQu'EiEil. 

The initial experimental and theoretical/numerical results on this system raised a number of very interesting 
questions: 

1) Is it possible to formulate a semiclassical theory of the tunneling current which involves only real (as opposed to 
complex) orbits within the well? If so, would such a theory only involve periodic orbits (and not all closed orbits as 
one finds in diamagnetic hydrogen), and would all or only a subset of the periodic orbits be important? 

2) The classical dynamics of electron motion within the well was observed numerically to have some unusual features. 
Certain regions of phase space became highly chaotic well before a KAM analysis would predict. Moreover, periodic 
orbits which connected the emitter and collector barriers were observed to disappear by backwards bifurcations, but 
it was rarely if ever possible to find the bifurcations in which these orbits were born. What was causing these unusual 
properties? 

3) Was the scarring of chaotic wavefunctions unusually common in this system? If so, what was the origin of this 
strong scarring, and why were only certain unstable periodic orbits reflected in the observed scarred wavefunctions? 

We attempt to give answers to all of these questions below, based on analysis of the classical dynamics and semi- 
classical methods. .Wp-note that this system has been extensively studied numericall^-both classically via the Suriape.- 
of-Section methodoull3, and quantum-mechanically by analysis rpf scaled spectrally and tunneling currentDi3'E2lEil. 
In addition to our own semiclassical treatment of the problems, recently there has been an independent work by 
Bogomolny and Rouben which presents two possible semiclassical tunneling formulas for this systerr£3. The first of 
these expressions neglects the structure of the emitter wassefunction and leads to a result which differs substantially 
from our periodic orbit forxnula, Eq. (^) below. In Ref£j the second was shown below to be equivalent to our Eq. 
(|l8|), first presented in Ref.fl, when similar assumptions were made about the emitter state (we present a similar proof 
and detailed comments below). ■— ■ n i— ■ 

Recent comparisons of this formula with exact quantum calculationaL-3 and with experimental datauo find good 
agreement. However, the periodic orbit formula was shown to lose accuracy and in some cases fail badly in narrow 
regionsEj near bifurcations, when rthe relevant periodic orbit does not yet exist or is not isolated. This failure is a 
manifestation of the "ghost" effectO. The relevance of this effect to the tilted well was first pointed out by Monteiro 
et al£3 who have recently shown that an adequate semiclassical description in this regime can be achieved by using 
complex dynamicsE-J Alternatively, we show below that one can use a representation of the tunneling current in 
terms of closed orbits (Eq. 14, section IIB), which involves only real orbits within the well and provides a remarkably 
precise quantitative description of the ghost effect (see section IVB). It is important to note that for the tilt angles of 
interest the spectrum cannot be calculated perturbatively from the 9 — limit, so that the only analytic option is a 
semiclassical approach, exploiting the fact that we are interested in states far above the ground state of the well, and 
the large level-broadening as discussed below. 



II. SEMICLASSICAL THEORY OF RESONANT TUNNELING 



A. Bardeen Approximation 

In standard, tunneling theory the current will be proportional to the transmission coefficient from the emitter to 
the collector^. This presents an immediate complication for a semiclassical formulation of the problem because 
the transmission coefficient involves the Green function across classically forbidden regions and hence would involve 
complex paths traversing the barriers. However another aspect of the physics suggests an alternative and simpler 
approach. The electrons are injected by tunneling from the emitter state into the well under high bias and far from 
equilibrium. Hence they will emit optic phonons and lose energy over a time scale of O.lps corresponding to 4-5 
collisions with the barriers, a time much shorter than the total tunneling time through the well which (from the 
measured tunneling current) corresponds to ~ 100 — 1000 collisions. Therefore there is a substantial level-broadening, 
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corresponding to several level- spacings within the well, and the tunneling is sequential, not coherent. Under such 
circumstances, in which correlations between tunneling into and out of the well can be ignored, one may employ the 
Bardeen tunneling formalismtH 

In Bardeen's "tunneling hamiltonian" approach, the system is treated as a sum of three isolated parts - the emitter, 
the quantum well, and the collector. Tunneling is considered as a perturbation, which introduces a nonzero coupling 
between these subsystems. The tunneling rate is then calculated using the Fermi Golden Rule, with the coupling 
matrix element expressed in terms of the wavefunctions of the isolated subsystems. One can then write a system 
of rate equations, expressing the balance of tunneling currents to and from the quantum well, and determine the 
nonlinear conductance of the system. 

If the transmission coefficients of the emitter and collector barriers are of the same order of magnitude, a substantial 
space charge is accumulated in the quantum well. As a result, the electron-electron interactions in the well become 
relevant, and it is not possible to connect the properties of resonant tunneling spectra to the classical mechanics of a 
single electron in the well. In order to probe the single-electron dynamics in the quantum well, the charge build-up 
in the well must be prevented. That can be achieved only when the width of the collector barrier is much thinner 
than the width of the emitter barrier, so that the tunneling rate from the well to the collector is much larger than the 
tunneling rate from the emitter to the well. In such a case we may set the occupation number of states in the well to 
zero in the rate equations for the current. 

In this regime the rate limiting step is the tunneling from the emitter to the well, and the tunneling current is 
therefore proportional to the corresponding tunneling rate (probability of tunneling per unit time). If there is only 
one energy level occupied in the emitter 2DEG (which is the case in the experiment at high magnetic fields B > 5T) , 
the tunneling current j is given by 



j = n e eW e 



(1) 



where n e is the surface concentration in the emitter layer. The Eq. ([!]) can be easily generalized to the case, when 
the emitter layer contains several occupied single-electron levels. 
The tunneling rate from the emitter to the well 



(2) 



where the summation is performed over all energy levels in the well, and the coupling matrix elementSlli is given by 
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with the wavefunctions ty e and corresponding respectively to isolated emitter and isolated well. 

The integration in (^) can be performed over any surface inside the barrier, including the boundary of the barrier 
itself. We choose to perform the integration over the "inner" boundary, which is classically accessible to the electrons 
trapped inside the well. The latter allows the semiclassical treatment of the wavefunctions of the isolated well using 
the Gutzwiller path integral approacbEj. 

The expression ([|) involves both the wavefunction of the isolated well ^ w and it's normal derivative S7^ w , taken 
at the surface of the barrier (z = 0) . In the limit of large height of the emitter barrier Uq corresponding to the 
experimental setupEa the normal derivative of the wavefunction does not depend on Uq, while for the wavefunction 

itself calculated at the plane of the barrier the dependence on Uo is crucial : ^ w (z = 0) ~ exp ^—y/Uo/£i \ , where 

the "injection energy" £j is defined as the difference between the total energy e and the potential drop across the well 
eV. Since both terms in Eq. ([}]) are of the same order of magnitude, Eq. ([?]) requires an adequate representation of 
both the well wavefunction and it's normal derivative. 

The standard semiclassical path-integral approach, which includes only real classical trajectories, does not distin- 
guish between infinite and finite potential barriers, and cannot be used for the calculation of the exponentially small, 
but non-zero value of the wavefunction at the plane of the barrier. This makes it complicated to apply directly the 
Eq. (||) to the semiclassical description of the resonant tunneling. Instead, we first reduce the expression (^) to the 
form, which involves only the normal derivatives of the wavefunction of the isolated quantum wellEl In the leading 
order in huj c /Uo, Si/Uo and eEa e /Uo (where a e is the width of the emitter barrier), the matrix element (||) can be 
expressed as 
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(4) 
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To obtain Eq. (Q) we neglect the effects of the electric field and the transverse component of the magnetic field in 
comparison with Uq in barrier region — a e < z < 0. Note, that for the untilted system (6 = 0) the approximation (^) 
gives exactly the same result as the standard WKB method. 

With the choice of the vector potential A = (—By cos(6>) + Bz sin(0))x the system possesses translational symmetry 
in theLjijdirection. Therefore, the dynamics within the well can be reduced to two degrees of freedom (see below) 
y, ZjloErtJ with an effective potential V (y, z). The well wavefunctions in Eq. (||) can be re-expressed in terms of the 
Green function of the isolated well, G(yi, zi = 0; y 2 , z 2 — 0; e): 



W e 



2ft 6 



d Vl dy 2 (ifc.O)*' (y 2 ,0) s 



d 2 G(y 1 ,z 1 ;y 2 ,z 2 ;e n ) 



dz\dz 2 



(5) 



zi=z 2 =0 



The emitter state in Ecu-.(jj) involves only the few lowest single-particle levels and can be calculated accurately 
using a variational approaches. The "injection function" ^> e (y, 0) is then a linear combination of the lowest few 



Landau level wavefunctions $ n , and is peaked at an injection point yi 
spatial extent of order the magnetic length, Ib = ^ft/eBcosO ~ ft 1 / 2 : 



(a e + I.13{ft 2 /m*eE) 1 / 3 )tan9, and has 
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(6) 



B. Semiclassical Approximation: The Closed Orbit Formula 



In the semiclassical approximation! 2 !], the derivative of Green function d 2 lZ2 G(yi, z\ = 0; y 2 , z 2 = 0; e) is determined 
by all classical trajectories connecting the points (yi, 0); (y 2 , 0) = (y — Ay/2,0); (y + Ay/2,0), and, as shown in 
Appendix A, can be expressed as: 
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=Z2=0 (2ttz) 3/2 ft 7/] 

where 5 7 (y — &ul2, 0; y + Ay/2, 0; e) = S 7 is the action of the classical trajectory0 indexed by 7, D 1 is the appro- 
priate amplitudeo, and pj and pj correspond to the initial and final momenta of the trajectory. Substituting (|7]) into 
(||), and introducing the effective broadening ft/T opt of the energy levels in the well, for the oscillatory part of W e -mi 
we obtain: 
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where t 7 is the classical propagation time of the trajectory 7, and the factor exp (— f 7 /T op t) represents the effect of 
phonon emission. This factor suppresses the coherent contributions of trajectories longer than r opt . 
We have also introduced in Eq. (||) the Wigner transform of the emitter wavefunction 



fa (y,p v ) = ft- 1 J dAy * e (y - Ay/2, 0) ¥* e (y + Ay/2, 0) exp (ip v Ay/ft) 



(9) 



which describes the distribution in transverse position and momentum of electrons injected into the well. Since ^ e (y) 
has a width ~ Ib, the integrand in Eq. (||) will be small for Ay > Ib ~ ft 1 / 2 . 

Assuming that M Landau levels contribute to the emitter state wavefunction, then the Wigner function 
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where the coefficients Uyk" can be obtained from the relation 
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where 



x L 1 . . , 2 . , 




(12) 



and is the generalized Laguerre polynomial. 

Consider first the integration over Ay in Eq. (|8|). Since Ay ~ ft 1 / 2 , in the semiclassical limit ?i — ► 0, one can expand 
the scaled action S~ / (y — Ay/2, 0; y + Ay/2, 0;e)/h in Ay, retaining only terms up to .second order, and perform exactly 
the resulting Gaussian integral. Alternatively, we can employ the approach of Berrju and perform this integration by 
stationary phase, which will initially lead to an expression in terms of non-closed orbits which satisfy the "mid-point 
rule" , and then re-express this answer to the same accuracy in H using Ay ~ y/K to arrive at the same result involving 
only families of closed orbits: 
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where M a = (mfj ) is the 2x2 monodromy matrixEJ, defined via the linearization of the Poincare map near the closed 
orbit a , and calculated at the contact point at the emitter barrier. 

The y-integration now involves the rapidly-varying phase exp[iS a (y, y)/h] for closed orbits beginning and ending 
at y. This phase is stationary for the periodic orbits, and we therefore expect the periodic orbits to determine the 
tunneling current. However, fw{y) varies on the same spatial scale ~ ft 1 / 2 , and we therefore cannot immediately 
perform the y integral by stationary phase (as is done to derive the trace formula for the total density of statestj). 
The failure of such an approach in this case reflects a more general principle of semiclassical theory; experimental 
observables typically do not depend only on the density of states, but also upon coupling matrix elements. The 
classical orbits which will contribute then depend on the nature and spatial extent of the coupling. In fact in the 
well-known problem of diamagnetic hydrogen, the "source function" analogous to fw(y) varies on a scale much smaller 
than the relevant phase, being effectively a delta- function, and the resulting semiclassical theory involves all closed 
orbits at the source pointcj and not periodic orbits. The current situation is intermediate between the absorption 
spectrum for diamagnetic hydrogen and a pure density of states measurement; a periodic orbit sum can describe the 
tunneling spectrum quantitatively, but a closed orbit formula (Eq. (|l4|) below) is more accurate near bifurcations. 
However the periodic orbit formula cannot be obtained by stationary phase integration at this pointLJ. 

Since the emitter state Wigner function fyy has the spatial extent of the order of the magnetic length I b ~ vft, and 
exponentially decays away from this interval, one can expand the scaled action S/fi in the powers of the deviation 
{y ~ Ui) ~ from the "injection" point y, (defined as the center of fw) and keep only the terms up to the second 
order. The integration over y will then yield: 
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where the summation is performed over closed orbits which start and end at the injection point y — y^. Generically 
there will only be one such trajectory in each topologically distinct family of closed orbits and this orbit will not 
be periodic, i.e. there will be a non-zero difference between the initial and fifftal momentum. In Eq. ( |i"4| ) = 
S a (yi,0. Hi, 0) is the action of this special closed orbit, n a is its topological indexed, and Ap a is the difference between 
the "initial" and "final" momenta: Ap^ — (p") y — {pf) y - In ( |l4| ) we have also introduced the following parameters: 
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where the coefficients Uk'k" , defined in (jl(]), represent the "transverse" dependence of the emitter state wavefunction, 
and are non- vanishing only for the few lowest Landau levels (k',k" ~ 1). The general expression for the matrix 
element Ck'k" is very complicated and will not given here. Here we only present the result for the most relevantraiililj 
families of the closed orbits, for which p" = 0: 
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where H k (x) is the Hermite polynomial, and £ a = (2 - Tr [M a ]) / (2 + Tr [M a ]). 

As noted above, typically there is only one trajectory which reaches yi out of each family of closed orbits and hence 
only one per family contributes to the tunneling rate in Eq. (|l4|). Let us call this orbit the injected orbit. Although 
the number of injected orbits increases rapidly with the increase of the orbit length, the contributions of long orbits 
with the traversal time t a longer than the relaxation time r opt are suppressed by the damping factor exp (— t a /T opt ). 
Since r op t generally corresponds to the time of a few traversals of the quantum well, it follows then that one may 
neglect all but the few shortest injected orbits and the sum in ( p^ ) is convenient for comparisons with experimental 
or numerical data. 

It is worth emphasizing again that the trajectories which contribute to the Closed Orbit Formula ([u]) are typically 
not periodic. The Closed Orbit Formula (COF) is therefore free from the problems which are encountered in semiclas- 
sical periodic orbit sums near bifurcations when quadratic expansions fail. We present comparisons of this formula 
with exact numerical results below in Section IV. B, after outlining our theory of the classical mechanics in the tilted 
well. These comparisons indicate that the COE-jjields an excellent quantitative description of the ghost effect, which 
is important near bifurcations in this systemt£fc3. The COF therefore represents a convenient alternative to the use 
of uniform approximationsEa and/or complex classical dynamicsEZI. The Closed Orbit Formula, Eq. ( |l4| ) is one of the 
main results of the current work. 



C. Periodic Orbit Formula 

As follows from (|lj), only injected orbits with the difference between initial and final momenta Ap y ~ \f%, can give 
substantial contributions to the tunneling current. Therefore for these important injected orbits, in the generic case 
(see Fig. || ) there exists one and only one periodic orbit (which has Ap y = 0) in a semiclassically small neighborhood 
near it. Using this proximity of each relevant injection orbit to a periodic orbit we can re-express the actions and 
momenta of the injection orbits in terms of the properties of their periodic neighbors (labeled by the index /z) as 
follows: 
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and for the tunneling current current we obtain 
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Here /i labels the periodic orbit, Sy = y — v r . 6p y — p y — (p^) y , the integer n M is the topological index@ of the periodic 
orbit, and the 2x2 monodromy matrixes M — (m^) is calculated at contact point at the emitter barrier. 

This relation between injected orbits and periodic orbits is generically valid because if Ap y is small for the injected 
orbit, then the function Ap y (y) for this family of closed orbits can be expanded linearly around yi and extrapolated 
to zero to find the nearby periodic orbit (see Fig. ||(a) ). However if that periodic orbit is in the neighborhood of 
a bifurcation this expansion will fail and Eq. (|l^) will become unreliable. The simple case of a tangent bifurcation 
is illustrated schematically in Fig. || (b). Here just before the bifurcation there are two periodic orbits which are 
close together, neither of which alone can adequately describe the contribution of the relevant closed orbits. Just 
after the tangent bifurcation there is no periodic orbit at all in the semiclassically small neighborhood of the injected 
orbit. There may however be another periodic orbit of the same topology far away from the injected orbit - see Fig. 
||(c) . In this case the linear expansion will yield a semiclassically large value for Ap y and it will not contribute to the 
sum, introducing no additional error, but missing higher-order contributions near the bifurcation. In such cases the 



periodic orbit formula (18) will lose accuracy but not fail completely as the resulting divergences are integrable (see 



Fig. Il2l below) . However for such cases the Closed Orbit Formula will remain quite precise and can be adapted for 
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the specific type of nearby bifurcation. In section IV. B we demonstrate this for the case of tangent bifurcations in 
which we calculate the exponentially small "ghost" contributions^] beyond the bifurcation. 

Eq. ( |l8| ) has the feature that it separates out the semiclassical summation over periodic orbits within the well from 

the contribution of the injection function fff (y,p y ) which cannot be calculated semiclassically and depends on the 
doping level in the emitter, as well as the magnetic field and tilt angle, thus it is convenient for a general quantitative 
description of the actual experiments. If the emitter state wavefunction has only zero-order Landau level component 
(C-k'k" ~ <5fc'o Sk"o), then the integration in ([l^) can be performed analytically. For the most relevant periodic orbits 
with p y = at the point of collision with the emitter barrier, we obtain: 
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where 
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h„ = n„ + 1 - Sign (Tr [M M ] + 2) . 



If the shift of the "injection point" yi from zero can be neglected, Eq. (|19| ) reducestll to the result, obtained by 
Bogomolny and RoubenLi 

The summation in ( |l8| ) is performed over all isolated periodic orbits, both stable and unstable. Near stable islands 
the motion is regular and we expect semiclassical quantization tq V i ic ld discrete energy levels (neglecting phonon 
broadening) and sequences of eigenfunctions localized on the island a 3 °n lc i. In contrast near unstable orbits the motion 
is chaotic and semiclassical theory does not yield discrete levelsEl This difference, can be displayed explicitly by 
performing exactly the summation over repetitions of the primitive periodic orbits in Eq. (flq), yielding 



Li' 

'cir 



dy I dp v fw(y,Py) 9e' (ViPy) 



(20) 



where A (tr, p) 



sinh(a") 



and the index + or — denotes stable or unstable orbits. The quantity h/r^g is an 



cosh (a") — cos(p) 

effective level-broadening which differs in the two cases. 

Note that the function A has a peak every time the semiclassical quantization condition S^ei) — 2-kTi{ti + n^/A) 
is satisfied, and these peaks become delta functions as t^ s — > oo. For stable orbits we find that = T opt , so if 
we neglect phonon scattering (r op t — > oo) we do recover perfectly discrete contributions to the tunneling current. 
In the stable case the argument of in Eq. (^p| ) is e\ = e — Tilu^ + (I + 1/2) which may be interpreted as the 
energy of longitudinal motion along the orbit. Due to the harmonic aprpoptipration the quantization of the transverse 
oscillations around the periodic orbit simply yields equally spaced levelal^'ES with spacing hco^ , where the frequency 
uj± = (fip/Tn) 4>fj, is the winding numberEEl and T M the period of the orbit. So the discrete energies at which tunneling 
occurs are the correct semiclassical energy levels of the well. 

The amplitude of each contribution is given by the coefficient functions gg in Eq. (|2p|), which are the Wigncr 
transforms of the harmonic oscillator wavefunctions corresponding to these transverse modes: 



9t 



(y, Py ) =(-!)% 2 



Q, 



exp 



(21) 



where L e is the Laguerre polynomial and = |2 + Tr [Af]| 1/2 |2 - Tr [M]\ lrz Q 



-1/2, 



Since the result (|2l| ) is based on the harmonic approximation within a stable island, we may only include modes 
up to £ m axj which is given by the ratio of the island area to ti. Phonon scattering smears out each of these discrete 
contributions to W osc over an energy range fi/T opt . 

In contrast, for unstable periodic orbits we find that t£ s = r op t/ (1 + ^A^Topt), where A M is the Lyapunov exponent 
near the orbit \i. Hence this time is finite and equal to 1/A^ in the absence of phonon scattering. Therefore, instability 
acts as a sort of intrinsic level-broadening, and each PO describes a contribution due to a cluster of levels. The peak 



of the function A in Eq. 
The weight functions 



GO) corresponding to the mean energy of the cluster is given by S^ie) = 2nh(n 



74). 



3e 



(y, Py ) = (-l)^ft Ue (2iQ^ exp (iQ^j (l + 



.sin(5 M /?i 



./2) 



sinh (T M /r^) 



(22) 
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are related to Wigner functions averagecBB over the eigenstates of the cluster. For each unstable PO the high-£ 
contributions are exponentially damped, and the main contribution to the tunneling rate is given by the £ = term. 

We now have a rigorous criterion for which periodic orbits contribute substantially to the tunneling current in Eq. 
(P0|). The injection function fw{y,Py) is centered on yi and p y — with widths ~ Ib,Ti/Ib in y,Py From Eqs. 
|2l| , |22| and the definition of we find that classical weight functions are are centered at y^, {py)^ with widths 

~ — y/2% \m^ 2 \ 1 4 — Tr 2 [Af M ]| 14 and h/l^ respectively. When the real and momentum space peaks of these 
two functions overlap the PO is semiclassically "accessible" and makes a substantial contribution to the tunneling 
current S3 



III. CLASSICAL DYNAMICS 



A. Model Hamiltonian 



The semiclassical formula for the tunneling current, Eq. (|20|), involves only classical periodic orbits within the well 
at energies below the top of the barriers. Therefore we can model the system by taking the barriers to be infinite 
perfectly-reflecting walls at z = 0, d. With our choice of gauge and coordinates we have 

H = {p x - eBycosje) + eBzsm(6)) 2 | P% | pj | ^ 
2m* 2m* 2m* 

+ U(-z) + U(z-d) (23) 

where the functions U ( U(z < 0) = 0, U(z > 0) = oo V,pepresent the infinite walls. Because H is independent 
of x, and p x can be .eliminated by a gauge transformation™], this problem can be simplified to a two-dimensional 
effective hamiltoniaro. The original hamiltonian involves four variable experimental parameters: B, E, 6, d and, as 
in the case of diamagnetic hydrogen, it is essential to express H in terms of the minimum number of dimensionless 
parameters. There are only two energy-independent classical length scales in the problem: the barrier spacing, d, 
and Id = (E/B)T C , where E/B is the drift velocity for perpendicular fields and T c is the cyclotron rotation period 
(T c = 2it/u} c ). Since the length-scale d will drop out of .the problem when the electron energy e < eEd = eV 
it is natural to scale all lengths by Z . When this is donea the effective hamiltonian is found to depend on three 
dimensionless parameters^ 2 ]: 6,(3 — 2vqB/E and 7 = e/eV, where vq is the velocity corresponding to the total 
injection energy e — m*v\j2 = + eV, where is the emitter state injection energy into the well. This energy is 
typically small compared to eV. Equally importantly, its ratio to the voltage drop is roughly constant in experimental 
conditions, implying 7 w 1.15 and constant. Hence for fixed tilt angle the scaled dynamics need only be studied as a 
function of one parameter, (3, which can be regarded as the scaled magnetic field. 



B. Poincare Map and the Origin of Chaos 



The most convenient surface of section to choose is the plane z = d corresponding to the collector barrier, as all 
trajectories within the well impinge on this plane (although in certain cases the emitter SOS may be useful as well). 
The most convenient phase-space coordinates for the system are the scaled velocities v x /vo,v v /vq (we note that for 
this system v x ~ y). When 6 = the (integrable) motion consists of cyclotron rotation perpendicular to the E field 
and uniform acceleration along it (perpendicular to the barriers) . Collisions reverse the ^-component of velocity and 
conserve separately longitudinal and cyclotron energy. When 6^0 the motion between collisions is still integrable, 
although now there is a drift velocity added to the cyclotron motion and the acceleration is along the magnetic field 
with magnitude ~ Ecosd. However collisions now exchange energy between the cyclotron + drift energy and the 
longitudinal kinetic energy. Since the tilt angle is in the y — z plane it is straightforward to see that there is no energy 
exchange at collisions when v y = and maximal energy exchange when v x — 0. This then is the origin of chaos in 
the system, the energy exchange between the two degrees of freedom depends up the phase of the cyclotron rotation 
at each collision, very much as in the Chirikov standard map (or kicked rotor )□, and just as in that case the system 
becomes chaotic as this frequency increases relative to the collision frequency (for non-zero tilt angle) - see Fig. ||. 
Hence as B increases or V decreases the system undergoes a transition to global chaos. A simple scaling analysis 
shows that for fixed.7 the classical parameter (3 ~ B/^/V, leading to parabolic curves of constant classical dynamics 
in the B — V planecl. n 

In order to derive a simple criterion for the transition to chaos in this system Shepelyansky and Stonea noted that, 
since 7 was close to unity, many trajectories would have insufficient longitudinal energy to reach to emitter barrier and 
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one could consider simply electron motion in the triangular well formed by the electric field and the collector barrier 
(the "single-barrier" model). In this approximation they showed that the dynamics reduces to an effective standard 
map with a smooth KAM transition to chaos and a chaos parameter given by K = 29(3. When K ~ 1 the last KAM 
curve of the system breaks and one has full exchange of energy between the degrees of freedom. However our recent 
work has shown that this approximation does not capture many crucial and interesting aspects of the physics, since 
those trajectories which do reach the emitter barrier are crucial, and including the second barrier changes qualitatively 
the dynamics in that region of phase space. The single-barrier model is a mathematically tractable and experimentally 
realizable KAM system; it is studied in detail in ref.B. Here we will only consider the double-barrier model relevant 
to current experiments and to our semiclassical tunneling formula above. 



C. Non-analytic Dynamics in the Double-Barrier Model 



The key new feature of the double-barrier model is that the presence of the second barrier introduces non-analytic 
behavior into the Poincare map. This map takes the values of v x ,v y at a collision and maps them to the values at 
the subsequent collisions. In the double-barrier model the map defines a closed curve in the SOS corresponding to 
those initial conditions which reach the second barrier with v z = exactly. The electron with the initial velocity 
inside this "critical boundary" reach the second barrier and get an additional "kick", and otherwise it will not. It 
is easy to show then that the map will have discontinuous derivatives on this boundary. As a result, just as in the 
stadium billiard or other systems described by non-analytic maps, the system in the vicinity of the critical boundary 
undergoes a sudden jump to chaos instead of a smooth KAM transition. However since trajectories sufficiently far 
from the critical boundary are confined away from the boundary by KAM curves, unlike the stadium one does not 
have global chaos for arbitrary tilt angle, instead one has a "chaotic halo" region surrounding the critical boundary 
(see Fig. |). This region becomes of substantial size well before the phase space becomes globally chaotic and hence 
has a significant impact on the transition to chaos. We will see below that all of the periodic orbits relevant to the 
semiclassical tunneling formula are born precisely at the critical boundary in the chaotic halo region. They are born 
in a novel type of bifurcation we call a "cusp bifurcation" which we will discuss in detail below. These bifurcations 
were very difficult to detect numerically from the SOS because they violate the bifurcation theorems for analytic 
hamiltonian maps. Understanding of the short periodic orbits and their bifurcations came from continuity arguments 
from the 9 = limit and from analytic results to be discussed below. 



D. Periodic Orbits at 9 = 



All periodic orbits can be easily classified at 9 = 0. First, there is a single isolate "traversing orbit" (TO) which 
has zero cyclotron energy and follows the line normal to the barriers, bouncing back and forth. Then there are an 
infinite number of families of "helical orbits" (HO's) with periods T given by the resonances between the cyclotron 
rotation with period T c and the longitudinal bouncing-ball motion with period Tl'. 

T = nT L = kT c . (24) 

In the collector Poincare map this would define a family of period-n orbits, all related to one another by rotation of 
the helix. An important point is that when 7 > 1 (the case we are interested in) there are can be two such families 
(or "rational tori") for each pair of integers (n, k). For simplicity assume that V and £0 are fixed and we are varying 
B. The resonance condition can then be rewritten as 

T L = (25) 
eBn 

where Tl is a function only of the fraction of the total energy which is in the longitudinal motion, el- The dependence 
of Tl on el is shown in Fig. |[ When El < eV then Tl increases as ^Jel\ however when el > eV then Tl begins 
to decrease because the electron reaches the emitter barrier and the bounce returns it faster to the collector barrier. 
Now consider, for any given resonance (n, k), increasing B from zero. Until B\ = it (m*E) 1 ^ 2 (2ed) k/n there are 
no solutions for the resonance condition, meaning that the cyclotron rotation is too slow for that particular resonance 
to occur (here we have used the explicit form of Tl{E,<1)). For B\ < B < B 2 = B\ (^7 + ^7 — l) there are two 
solutions corresponding to two different families of helical orbits with the same Tl; one of which has longitudinal 
energy greater than eV and reaches the emitter barrier, the other of which has El < eV and doesn't reach the second 
barrier (see Fig. ||). We will refer to these two types henceforth as "emitter" and "collector" orbits. As B is increased 
further the emitter orbit must give up cyclotron energy to remain in resonance, while the collector orbit must gain 







cyclotron energy in order to remain in resonance. Eventually at B = B2 the emitter orbit has zero cyclotron energy 
and is degenerate with the traversing orbit, which exists at all magnetic fields for each energy. For B > B2 the emitter 
family of periodic orbits ceases to exist, but the collector family persists for all B > B-y. 

In the surface of section this evolution look as follows (Fig. g): at B\ a circle of degenerate period- n orbits appears 
exactly at the critical boundary (i.e. reaching the emitter barrier with v z — 0. As B increases this splits into two 
circles, one with radius greater than the critical boundary (the collector torus) and one with radius smaller (the emitter 
torus). As B tends to B\ the radius of the inner circle shrinks to zero and the torus disappears. The conclusion is 
that at 9 — all periodic emitter orbits (except the traversing orbit) exist only for a finite interval of magnetic field 
(and a similar conclusion holds if V is varied as well). By continuity the same must be true for 9 =/= 0. 

E. Effect of Non-zero Tilt on Periodic Orbits 

Tilting the magnetic field makes the dynamics non-integrable and as usual all families of periodic orbits are de- 
stroyed, being replaced by an equal number of stable and unstable isolated fixed points with the same periodicity in 
the map. Isolated periodic orbits (in this case thg traversing orbit) survive the perturbation, and the resonances of 
these orbits transform into its bifurcations. In ref.tl we have given an exhaustive analysis of the classical periodic orbit 
theory in this system. Here we will simply consider briefly the case of period-two orbits which is both illustrative 
of the general situation and most relevant to the experiments. In particular we will consider the period-two orbits, 
which occur at the lowest values of magnetic field (these orbits at 9 = make 1/2 a cyclotron rotation per collision 
with the emitter). There are four such orbits, with different topology, the (0, 2)a,(2, 2)a, and (0,2)y, (2,2)y orbits 
(see Fig. 0). The first and second numbers in this notation represent the numbers of bounces per period with the 
emitter and with the collector respectively, and the subscript describes the "topology" of the orbit. 

Since the emitter orbits have different topology they cannot be created together e.g. the (2,2)a must be created 
with the (0, 2)a and the (2, 2)y with the (0, 2)y. Since (0, 2)a,(0, 2)y are collector orbits and (2, 2)a,(2, 2)y are emitter 
orbits, this can only happen precisely at the critical boundary (just as for the unperturbed case); the bifurcation at 
which they are born has the non-analytic properties mentioned above which lead to a cusp at the bifurcation point 
in the bifurcation diagram as shown in Fig. |g. In correspondence to the evolution of tori at 9 = shown in Fig. [|, 
now as P (or magnetic field) is increased the emitter orbits fixed points move inward in the SOS (but now at slightly 
different rates for the two different emitter orbits), eventually annihilating with the period-one traversing orbit in 
successive backwards period-doubling bifurcations. This scenario is general for all emitter orbits: birth at the_critical 
boundary in a cusp bifurcation; death in a backwards bifurcation in the neighborhood of the traversing orbiio. 

There is an important sub-plot in this scenario. Once the symmetry of the problem has been broken by non-zero 
tilt angle it is not possible for a (2, 2)y and (0, 2)v orbit to be born together, because the two "arms" of the V orbit 
can no longer be equal in length. Hence the two (1,2) orbits are born first, one of which has an arm reaching the 
emitter barrier and one of which does not. Withjncreasing fj the emitter orbit evolves into a (2, 2)y orbit, typically 
not directly, but through a bifurcation "cascade"Ej of the type shown in Fig. ||. In such a cascade a (1, 2) orbit born 
in a cusp bifurcation with a (0,2) will disappear in a tangent bifurcation with another (1,2) orbit born in a cusp 
bifurcation with a (2,2). Thus the first (1,2) orbit is not directly connected by bifurcation to the (2,2), but only 
indirectly as shown in Fig. |8|. 

Such bifurcation cascades connected to a resonance of the traversing orbit are a generic feature of this system, for 
period- N orbits. The only difference is that for the orbits with many bounces in the corresponding bifurcation cascade 
the "connection" of the (0, N) and (N, N) orbits, required by continuity to the behavior at zero tilt angle, is a "multi- 
step process", which includes all the "intermediate" periodic orbits of the sequence (1, N), . . ., (N — l,N). Thus, 
despite the increasing complexity, all bifurcation cascades in this system can be analyzed within a single framework. 

Each of these bifurcation cascades terminates finally in a resonance with the periodi-one traversing orbit. The 
analytic expressions for the period and stability of all period-one orbits, derived in Ref.til, provide a straightforward 
description of each of these resonances, and therefore allow a direct and unambiguous identification of all of the 
primary bifurcation cascades. 

F. Cusp Bifurcations and Stability Analysis 

The existence of both emitter and collector period-one helical orbits at 9 = 0, which appear together at the critical 
boundary, implies that for 9 7^ their counterparts must be produced in cusp bifurcations. The analytic analysis 
of the period-one orbits allowed us to obtain a detailed description of these cusp bifurcations. For example we were 
able to show that these orbits do not tend to marginal stability at the cusp bifurcation, in contrast to the behavior 
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required of analytic maps. This led to some more general results on the properties of all CB's, which are as follows: 1) 
all cusp bifurcations appear at the critical boundary, where the Poincare map is non-analytic; 2) any cusp bifurcation 
produces two periodic orbits, which differ by one "point" of contact with the emitter wall; depending on the topology 
of these orbits, they differ either by one or by two bounces with the emitter barrier per period; 3) out of these two new 
orbits, the one which has more bounces with the emitter wall (the "more connected" orbit), is "infinitely unstable" 
(Tr[M] — * oo), while it's partner might either be stable or unstable. 

The calculation of the trace of the monodromy matrix for the V-type period-2 orbits, presented in the Fig. |j] , 
illustrates the latter point. Note the divergence of the Tr[M] for the orbits (1,2) + (the "more connected" partner 
of (0, 2)v) and (2, 2)y (the "more connected" partner of (1,2)~). Meanwhile, their partners appear either stable 
((1,2)-) or unstable ((0,2) y ). 

We hypothesize that the presence of cusp bifurcations with these properties is a general feature of hamiltonian 
maps with hard wall potentials when the orbits can either hit or miss one part of thc_hard wall. Therefore the 
same phenomenology should exist in certain billiard systems (e.g 1 — the "annular billiard"E3), or in billiards in strong 
magnetic ficldE3, the possibilities we are currently investigating^. This phenomenology has importance for the 
quantum properties of the system, because it leads to special 'metastable orbits" as we now discuss in the case of the 
period-two orbits in the tilted well. 

As the effective chaos parameter /3 increases, the collector orbit (0, 2)y becomes increasingly more unstable, as 
expected for a periodic orbit of a system undergoing the transition from integrability to chaos. However, the other 
"least connected" orbit of this bifurcation cascade, the orbit (1, 2) _ , shows a dramatically different behavior. Initially, 
with the increase of (3 it loses stability, as expected. However, this orbit will later have to annihilate with the orbit 
(1,2) + in a backwards tangent bifurcation. Since the tangent bifurcation can only involve two marginally stable 
orbits, the initially stable (1,2)~ can not develop strong unstability and must stay either stable or only marginally 
unstable in it's whole interval of existence. 

Note, that this continuity argument applies to any "intermediate" period orbit, which appears as the "less con- 
nected" partner in a cusp bifurcation, and annihilates in a tangent bifurcation. Therefore, there must be at least one 
such metastable orbit in any primary bifurcation cascade of this system. This "metastability" of a certain subset of 
the classical periodic orbits in the quantum well in a tilted magnetic field has profound consequences in its quantum 
spectrum. 



G. The Origin of Strong Scarring 

It is now well known that certain quantum wavefunctions of classically chaotic systems can strongly localize near 
unstable periodic orbits. Such wavefunctions_arje-|termed "scars" and their properties have been extensively studied 
and well understood during the last decadesE3~El. The scarred wavefunction are most likely to be associated with 
the periodic orbits which are not too unstable. A sufficient condition of the strong scarring is AT % 1, where A is 
the instability (Lyapunov) exponent associated with the orbit, and T is it's period. Typically, in chaotic systems, the 
periodic orbits quickly become increasingly more unstable, as the classical parameter driving the system to chaos, is 
increased. Such orbits can strongly scar the wavefunctions only in a small interval of classical parameter space in 
which they are not too far from stability. 

Therefore, it came as a surprise, when a persistent scarring by,-the same orbit in a large interv&Lof applied voltage, 
was observed both numericallyQ and (indirectly) experimentallyll Although it has been arguedc3 (with some merit 
in our view) that the experimental results are in a regime of low quantum numbers for which the concept of scarring 
becomes ambiguous, extensive quantum numerical calculations confirm a strong periodic scarring at high quantum 
numbers as well. This phenomenon is naturally explained by the "metastability" of the "intermediate" periodic orbits 
in this system. Since these orbits remain close to stability over their whole interval of existence, they produce strong 
scars in the quantum wavefunctions in a large interval of the applied voltage or magnetic field. This behavior is 
illustrated in Fig. |l^, which represents the calculation of the scar "strength" (defined as the value of the Husimi 
projection of the quantum wavefunction, taken at the fixed point the periodic orbited) due to the (1, 2)~ orbit. 

As we already pointed out, there is at least one "metastable" orbit in each of the primary bifurcation cascades. 
The semiclassical theory therefore implies, that each of these orbits should also produce strong scars in quantum 
wavefunctions. Indeed, a detailed analysisu of the wavefunctions in the quantum well uncovers strong and persistent 
scarring due to the "metastable" period-3 orbit (1,3) and period-5 orbit (1,5). 
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H. "Exchange Interaction" between Bifurcation Cascades 



An interesting feature of the classical dynamics in the quantum well in a tilted magnetic field, noted in RefJl is the 
exchange bifurcations between the orbits of the different primary bifurcation cascades. Fig. [ll] depicts the exchange 
bifurcation between y-type periodic orbits of two different period-2 cascades. For convenience of the further analysis, 
we represent the relevant orbits by the x-components of the velocity v x oc y at the emitter barrier (by definition, such 
diagrams do not show the collector orbits). As the tilt angle approaches the critical value of 9* w 30°, the orbits 
(1,2)~ and (1,2) + of the first cascade, and the orbit (1,2)+ of the second cascade come close to each other, and at 
the critical angle exchange partners, so that for 9 > 9* the orbit (1, 2)+ is now paired with (1, 2)~ and belongs to the 
first cascade, whereas the orbit (1,2) + is paired with (1,2)~ and belongs to the second cascade. As the tilt angle is 
increased further, the two reconfigured bifurcation cascades move away from each other - see Fig. [lip . 



IV. SEMICLASSICAL TUNNELING SPECTRA AND COMPARISON WITH EXPERIMENT AND 

NUMERICAL CALCULATIONS 



A. Periodic Orbit Theory 



Since the motion of the fixed points of the relevant periodic orbits near the exchange bifurcation strongly affects 
their accessibility for tunneling, the exchange bifurcation dramatically manifests itself in the tunneling spectra. This 
is demonstrated in Fig. [fl|, where the gray-scale regions (calculated at B = 8T) represent the semiclassical width Z M 
of these wavefunctions localized near the corresponding periodic orbits. The width of the injection function fw(y) 
is represented by the region between dashed lines. Whenever these regions overlap for some value of (3 the orbit is 
accessible and Eq. ( |20| ) predicts that a peak-doubling-cegion will appear in the B — V parameter space along the 
parabola corresponding to that value of (3 (see Fig. |l2|)c-3. 

At 8 — 29° the accessibility intervals for the (1, 2)~ orbit and the (1, 2)+ orbit cover the entire interval 4.3 < (3 < 10.9 
overlapping briefly around (3 = 7.5 n Thus one expects a large region of peak-doubling in the B — V plane with no 
gaps as observed in the experimentn. Since the periods of the relevant orbits are almost degenerate, these distinct 
periodic orbits will produce the oscillations of the tunneling current of almost the same frequency. Note, that the 
orbits (1, 2)~ and (1, 2) + of the first cascade are slightly better coupled with the emitter state, that the orbit (1, 2)+. 
Therefore, in the I-V, we expect the amplitude of oscillations to be higher in the high- voltage region (/? ~ 1/y/V). 

As 9 = 31°, immediately after the exchange bifurcation, the relevant orbits of the newly- formed "first cascade" 
start to move away from the accessibility region. As a result, the amplitude of the oscillations of the tunneling current 
is now higher in the low-voltage region, and we expect an abrupt amplitude change around [3 w 7.3, (V ~ 0.44V for 
B = 8T). This is observed clearly in the experimental data of Fig. |l2|(a),(d), where the amplitude change is compared 
to theory. The theory is found to predict the ratio of the amplitudes between the low and high voltage regions (above 
0.45V) with only 25% mean error. 

At s» 34°, the periodic orbits of the first cascade are no longer accessible. At the same time the (1,2) + orbit is 
still highly accessible at low voltages (high (3). Thus we expect the high voltage oscillations to disappear at 9 = 34° 
while the low voltage oscillations persist. This is seen clearly in Fig. [l^, again in good qualitative and quantitative 
agreement with theory. _ 

In Fig. [l2[ we plot the peak position data of RefE in the entire B — V plane against the semiclassical prediction based 
on the contributions of these periodiCpOrbits. Again, the comparison shows a good agreement between semiclassical 
calculation and the experimental datac3. 



B. Beyond Periodic Orbits: The Closed Orbit Theory 

As was demostrated in the previous section, the periodic theory can be successfully used for the quantitative 
description of the tunneling spectra over a wide range of system parameters. However, the periodic orbit approach 
may fail near the bifurcations, when the relevant orbits are not isolated, and the second-order expansion, used to 



obtain the POF, is not justified. A manifestation of this can be seen in Fig. 12i,d, where one sees a spurious dip 
in the semiclassical trace near V = 0.3V which has no analog in the experimental I — V curve. This dip is precisely 
due to the occurrence of a bifurcation near this voltage. Also, the POF provides no description for the "ghost effect" 
near tangent bifurcations, when the contribution of the particular orbits which disappear in the bifurcation remains 
in the spectrum even after their disappearance^. 
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Generally, the regions where the POF is inadequate, are semiclassically small, i.e. the corresponding intervals 
AV, AB ~ h. And in the case of the ghost effect, the ghost contributions decay exponentially away from the 
bifurcation point. However, due to a relatively large value of the effective h in the experimenturoQ (particularly in 
Rcf.El), the measured turuaeling spectra do show clear features related to these effects. In particular, as it has been 
shown by Monteiro et al.td that a substantial part of the peak-doubling region at 8 = 27° is due to the ghost of the 
(1, 2) orbit, and the period-1 oscillations at high magnetic field at 9 — 11° are related to the ghost of (1, 1) orbit. 

Previously the "ghost" contibutions have been calculated using an extension of the actual classical dynamics to the 
complex space. In this approach, after a real periodic orbit "dies" in a bifurcation, a counterpart orbit continues to 
exist in complex space-time. The imaginary part of the action of this complex periodic orbit describes the exponential 
damping of it's contribution in the density of stateaij. 

Very recently, it has been shown, that complex dynamics can be successfully used for the quantitative description 
of the tunneling current)!! A full stationary phase condition, when applied directly to Eq. (|^), leads to non-periodic 
and complex orbits. These trajectories were shown to smoothly "interpolate" between successive relevant periodic 
orbits, also providing an excellent-quantitative description of the tunneling current in the "ghost" regime. 

However, the approach of RefJiZI does not possess one of the important advantages of the semiclassical theory 
described in the pr£sent work - the separation of the" isolated" subsystems of the quantum well and the emitter state. 
In the approach oOa different emitter state leads to a completely different set of complex trajectories, which enter 
the final expression for the tunneling current. 

The Closed Orbit Formula derived in the present paper ( Eq. (14) ) has the advantages that it only involves real 
orbits within the well, which depend on the emitter state only weakly (through the location of the injection point 
Hi). It is therefore of interest to see if a semiclassical formula with this rather different structure can also describe the 
ghost effect in the tunneling spectra quantitatively. 

The relations (|l3|) , (fl4}) provide a complete description of all the ghost contributions in the tunneling spectra. For 
simplicity here we will apply them to self-retracing orbits (the most relevant to experiment) and introduce a further 
approximation to re-express the injected orbit in terms of the closed orbit with the minimum Ap y . 

As it was argued in Section II. C, an injected orbit can substantially contribute to the tunneling current only if 
the change of momentum Ap y after returning to the injection point is small. Consider such an injected orbit when 
it is near a tangent bifurcation as shown in Fig. p[ Just before the bifurcation there exist no periodic orbit in 
a semiclassically small neighborhood of the injected orbit, but there does exist a unique nearby orbit of minimal 
Ap y . This orbit is determined by the classical dynamics in the well and not by the injection point, so we prefer to 
express the semiclassical tunneling current in terms of its properties, rather than those of the injected orbit. Let its 
point of contact with the emittee-.be at yo, where \yo — yt\ is semiclassically small by assumption; than the change of 
momentum of the injected orbitCJ can be expressed as 



A Py = 2 Py + K ( Vi - yo) 



Kit 



2 Py 



dy 2 



(26) 



Assuming that the electrons in the emitter (or collector) are in the lowest Landau level, for the corresponding ghost 
contribution to the tunneling rate we obtain: 



Wgho 



A 



where 



gh ex P 



to \ 

Topt J 



So 



(27) 



^ c h = A exp (-C 2 



(28) 



2/7 



1/4 



exp 



{PyiB/ny 



(29) 



(py + Koiy.-yof) 
1 + ( Kq {yi - y ) 



1/2 



(30) 



m®j , to and St) are respectively the monodromy matrix, the traversal time and the action of the closed orbit with 
minimum change of momentum Ap y , and the </>o is the corresponding phase, which, in addition to the constant part, 
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related to the topological index of the orbit, has also a smoothly varying component. Although it is straightforward 
to obtain the latter from (|l4|), the result is somewhat complicated, and is not presented here. Note, that the classical 
input for the E q. ( j 27| ) depends again only on the injection point and not on the particular emitter state. 



Thus in Eq. (|27|) we have derived a relatively simple explicit formula for the ghost contributions in the tunneling 



spectrum in terms of the classical properties of the closed orbit with minimal values of Ap (= 2p y for a self-retracing 
closed orbit). Note that the strength of the ghost contribution is exponentially small in (Ap/ps) 2 , whereps = U/Ib ~ 



Now we proceed with a detailed numerical test of the "ghost formula" (|27|). Note that in principle the expression 
( p7| ) can describe the tunneling rate both from the ground Landau level in the emitter to the well or from the well 
to the ground Landau level in the collector. Although the actual experimental tunneling current in the incoherent 
regime does not depend on the latter quantity, either can be used to test the validity of Eq. ( p7| ) in the ghost regime. 
In fact the tunneling rate into the collector is a more convenient object for numerical computation since it is non-zero 
for energies less than eV and thus can be expressed in terms of lower-lying eigenstates for interesting value of eV. 
Effectively then we will test the formula in the single-barrier model which we studied in detail in referenced for which 
the classical dynamics depends only on (3. The first (as a function of (3) relevant bifurcation in that model is the 
tangent bifurcation of the period-one orbits which we designated (0, 1) + W and (0, (for a detailed description 
of the corresponding classical dynamics, see Ref.0, Section III.C). We will now consider the "ghost" related to that 
tangent bifurcation at 9 = 11°. Although for any given system parameters yi and yo are fixed in Eq. (p7|), for a 
complete test of the formula we treat yi as an independent variable and study the amplitude A s ^ h as a function of 
both yi and (3 (~ ^/e). 

In Fig. H we show the comparison of the exact and semiclassical results for Aq . With no fitting parameters the 
semiclassical calculation (represented by a dashed line in Fig. |l^) is in a good agreement with exact data (open 
circles). However, essentially perfect agreement can be obtained when one introduces a single fitting parameter by 
means of the definition ^4q c (el) where the "longitudinal energy" el is defined as El = e — kTiuj Ci with the K ^ 1 being 
the fitting parameter. 

The quantity k,Tiu c can be interpreted as the energy of the motion transverse to the direction of the periodic orbit. 
As was shown in Section II. C, a careful consideration of the higher-period orbits, related to the repetitions of the short 
trajectories, leads to the replacement of the total energy by it's longitudinal part in the phase of the corresponding 
contribution to the tunneling rate. It is natural to assume, that such a replacement should take place in both the 
phase and the amplitude. However, because the transverse energy e± ~ Ti, the correction to the amplitude is of the 
next order in H, and is not captured by the semiclassical theory. The calculation, presented in Fig. [l3] gives a clear 
indication that the amplitude should also, in fact, depend on the longitudinal energy. 

As follows from (^7|), neglecting the slowly varying prefactor, the amplitude depends only on the "scaled deviation" 
from it's maximum value £: 



In Fig. [l4| we plot the scaled amplitude as a function of £ for eight different values of (3, and a good agreement is 
found with the predicted form (|jl|). Note, that since the formula ( f27| ) is, in turn, also based on the expansion near 
the contact point yo of the "minimum Ap y closed orbit" , the deviation of the exact data from (|3l]) should increase 
for larger £, again in agreement with the calculation presented in Fig. 



In the work which we have concluded with the present paper, we have developed a complete semiclassical theory of 
resonant tunneling in quantum wells in a tilted magnetic field. The theory was shown to be in good qualitative and 
quantitative agreement with the experiments of Muller et alo, and explains many puzzling features of the tunneling 
spectra. In particular we have explained the phenomenon of persistent scarring of wavefunctions in terms of the 
properties of cusp bifurcations and bifurcations cascades, which lead to metastable periodic orbits. We expect these 
concepts to extend to other systems and for the same persistent scarring to be found in such systems in the appropriate 
parameter ranges. Furthermore we have demonstrated explicitly the possibility of describing quantitatively the ghost 
contributions of vanished periodic orbits in terms of a semiclassical theory of real closed orbits. 

The authors wish to thank G. Boebinger, D. Shepelyansky, T. S. Monteiro, T. M. Fromhold, H. Mathur, G. Hack- 
enbroich and J. Delos for helpful discussions. We further thank T. S. Monteiro for communicating the results of recent 
work on the semiclassical formulas and their extension to complex orbits prior to publication. This work was partially 
supported by NSF Grant No. DMR-9215065. 




Vh. 



A™ h (yi,P)~exp -£(y 4 ,/3) 2 



(31) 



V. CONCLUSIONS 
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APPENDIX A: SEMICLASSICAL GREEN FUNCTION NEAR A RECTANGULAR BARRIER 



The equation (^), used for the calcluation of the tunneling rate from the emitter to the well, requires the knowledge 
of the normal derivative of the Green function a~ , Z2 G (j/i, z\ = 0; y 2 , z 2 = 0). Because of the apparent discontinuity 
of the potential exactly where the Green function is evaluated, the semiclassical approximation for this object should 
be performed with care. 

First, we consider the first derivative d zl G(y\,Z\ = 0;y 2 ,z 2 ), which we define as 



9G(yi,zi=0;y2,z 2 ) ,. dG (yi, z x ; y 2 , z 2 ) - G(yi,0; y 2 ,z 2 ) 

a = llm 

az\ zi^o z\ 



(Al) 



The value z\ — > +0 should be classically small (e.g. Z\ ~ U), so that the change of the smooth part of the total 
potential V(y, z) at the distance z\ is small compared to the kinetic energy of the electron. .— . 

When z\ ^ 0, the semiclassical Green function can be expressed as a sum over classical trajectoriesEl (indexed by 
v) connecting the points (yi, Zx) and (y 2 , z 2 ): 



G 



2tt 



{2i:iTi) 



3/2 



exp 



(A2) 



where D is the (semi) classical amplitude, S is the action integral and /i is the topological index of the trajectory v. 

When the starting point [yi,Z\) is near the emitter barrier, the Green function can be expressed in terms of the 
classical trajectories starting at the emitter barrier - i.e. from the point (y±, 0). For each such trajectory j there are 
two trajectories j% and j 2 [starting from (yi, 1)] with different number of collisions with the barrier, which collapse to 
j when z\ — > - see Fig. |l^. In the semiclassical limit we have: 



S h = Sj -p^zx, S n = Sj +p i z 1) z 1 , 



(A3) 



where p^ 1 ^ the electron momentum at the starting point, and the extra 2 in the expression for fij 2 corresponds the 
phase increment of 7r due to a n additional (with respect to the trajectory j) reflection from the hard wall. 
Substituting (A3) into (A2), we obtain: 



G(yi,zi;y 2 ,z 2 ) 



{2-KiK) 



3/2 



\Dj\ sin 



' (i) 1 

Pz Z\ 



exp ( - Sj 



2 J 



(A4) 



and therefore the first normal derivative 
dG (yi,zi;y 2 ,z 2 ) 



dz\ 



01=0 



(2«) 3/2 a 5/2 , 



2 J 



(A5) 



Using a similar approach, one can show, that the second derivative 
d 2 G(y 1 ,z 1 ;y 2 ,z 2 



dz\dz 2 



!1=Z2=0 (2iri) 3/2 h 7/ '- 



■E^ 1)r,(2) 



Pz 



Dj exp - 5. 



infij \ 



(A6) 



where pi and pf are respectively the initial and the final momenta. 
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FIG. 1. Schematic of the geometry and coordinate system used in this work 




y, arb. units 



y, arb. units 



y, arb. units 



FIG. 2. The change of momentum Ap y for the closed orbits of a particular family, as a function of "starting coordinate" y 
The dashed lines show the linear approximation for Ap y , used in (^), ([t7|). The arrow indicates the "injection point" y = y t 
If at y = yi the value of Ap y is small, the linear approximation is justified (a), but fails in the vicinity of a bifurcation (b),(c) 
In (b) the dotted and solid lines respectively represent the behavior before and after a tangent bifurcation. 
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FIG. 3. The evolution of two classical trajectories, started from the same point at the emitter barrier with almost identical 
velocities, for (a) ui c t < 1, and (b) uj c t ^> 1, where r is the well traversal time. The black arrows represent the velocities 
immediately after collision with the collector barrier; the inset shows the directions of the magnetic and electric fields, and the 
coordinate system. Note the large amplification of small differences when the cyclotron rotation rate is large compared to the 
traversal time. This is the origin of chaos in the system, as discussed in the text. 
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FIG. 4. The "chaotic halo" near the critical boundary (thick solid line), created by nonanaliticity of the Poincare map, seen 
at the collector barrier surface of section at (3 = 2, 7 = 1.17 and 6 = 30°. Note that the halo region in which trajectories can 
hit the emitter barrier intermittently is highly chaotic while the inner and outer regions, where trajectories either always hit or 
always miss the emitter, are still quite regular. All the periodic orbits relevant to the tunneling current are born at the critical 
boundary and then move inwards. 
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B < B 




eV/£ 1 8 L /8 

FIG. 5. The dependence of the traversal time of longitudinal motion tL on the fraction of the total energy which is 
longitudinal, £l/e. The dashed lines represent the three different regimes discussed in the text, B < Bi, Bi < B < B2, B > B2 
and the corresponding borders B = Bi and B = B%. The insets on the right schematically represent the different periodic 
orbits in each regime. 
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FIG. 7. Examples of the different types of period-2 orbits in y-z projection. At higher magnetic fields similar orbits arise 
with one or more extra turns around the magnetic field direction. 
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FIG. 8. The bifurcation diagrams of the period-2 A-orbits (a) [6 = 11°] and F-orbits (b) [9 = 17°]. Each of the period-2 
orbits is represented by it's velocity [v x for A-orbits and v y for V-orbits] at the collisions with the collector barrier, with the 
velocity at the first collision shown by solid line, and the velocity at the second bounce shown by dotted line. The dashed line 
represents the period- 1 orbit (1,1). 
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FIG. 9. The trace of the monodromy matrix for different period-2 orbits at (a) 8 — 17° and (b) 9 — 28°. All these 
period-two orbits appear in cusp bifurcation. One cusp bifurcation produces the A-orbits (0,2)a and (2, 2)a, another cusp 
bifurcation generates the orbits (1, 2)~ and (2, 2)v, and yet another cusp bifurcation produces (0, 2)v and (1,2) + . Note the 
discontinuity of the trace of the monodromy matrix of the "more connected" partners (2, 2)a, (2,2)y and (1,2) + . The single 
"metastable" orbit (1,2)~ leads to a strong scarring, it's interval of existence substantially increases with the tilt angle. 
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FIG. 10. A wavefunctions scarred by unstable (1,2)~ orbit with y-z projections of orbit superimposed (a), and the "Scar 
strength" H vs. the scaled action S(e n )/h of the (1,2)~ orbit which scars the eigenstate of energy e n . The arrow indicates 
the value of f3 for the tangent bifurcation, which give birth to the periodic orbit. Note that when increasing action (or 
energy) the tangent bifurcation at higher /3 ~ l/v^ | e fPf rs a * the lower action side. The peaks of the scar strength below 
the tangent bifurcation are due to the "ghost effecf'LJ'Lj. Scaled action below the bifurcation points is obtained by linear 
extrapolation of the (approximately linear) function S(e)/h. "Scar strength" H is defined as the value of the Husimi function 
H(y°,Py) = J dy dp y W„(y,p y ) exp (— (y — y°)/l% ~ ^%(Py ~ P») 2 /^) calculated for the normal derivative of the wavefunction ( 
W n (y,p y ) — h~ x J dAy d z ^ (y — Ay/2, 0) d z ^ (y + Ay/2, 0) exp (ip y Ay/h) ) at the location of the fixed point (y° oc v x ,py oc v y ) 
of the unstable periodic orbit. 
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j8 (a) j8 (b) jg (c) 

FIG. 11. Bifurcation diagrams for the four period-two orbits relevant to the peak-doubling in the interval 8 = 29° — 34° 
(see text). Thick lines indicate y-coordinate of the single collision with the emitter; these lines coincide at bifurcations. 9 = 29° 
(a) indicates behavior before the exchange bifurcation; (9 — 31° (b) and 34° (c) after. Shading represents the localization 
lengths Ifj, associated with the relevant orbits [(1, 2)~ and (1,2)+ for (a), and (1, 2)~ and (1,2) + for (b),(c)]. Region between 
the dashed lines denotes semiclassical width of emitter state; overlap indicates a large contribution to the tunneling current. 
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FIG. 12. Figs. (a),(d): resonant tunneling I — V traces for 9 = 31°, # = 34° at B = 8T. Trace (1) is raw experimental 
data, trace (2) is same data, filtered to retain only period-two oscillations, trace (3) is semiclassical prediction from Eq. (po[). 
The modest discrepancies in the shape of the envelope of the amplitude of the oscillations are due to the inaccuracy of the 
quadratic semiclassical theory near the bifurcation which occurs around 0.3V at 8T. Figs. (b),(e): peak positions vs. voltage 
and magnetic field, determined from multiple sets of experimental I — V data at 6 — 31°, 34°. Figs. (c),(f): semiclassical I-V 
oscillations for same, note gray-scale indicates relative amplitudes, not just peak positions. Note disappearance of high vo ltag e 
oscillations at 8 — 34° due to movement of (1, 2lj~ orbit away from accessibility after the exchange bifurcation (see Figs. |ll|b, 
|rf]c and text) [The Figure is reprinted from Ref.u]. 
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FIG. 13. The dependence of the amplitude of "ghost" oscillations on /3 for 9 = 11°. The open circles represent the results 
of the exact calculation. The dashed line corresponds to the semiclassical formula (p9[). The solid line corresponds to the 
semiclassical calculation, which takes into account the transverse energy kTiuj c , with k = 0.55. This is in a good agreement with 
a simple estimate of k « 0.42 from the transverse energy of the stable (0, 1) + orbit, which is soon to be born in the tangent 
bifurcation. 
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FIG. 14. The dependence of the scaled amplitude on the scaled deviation £ from it's maximum value. The solid line 
represents the functional form (|3l|). Different simbols represent the results of the exact calculation for different /3's, with the 
corresponding values shown in the inset. Since the result (^l|) is based on the expansion near £ = 0, the deviations from the 
semiclassical result increase for larger £. 




FIG. 15. The classical trajectories, which contribute to G (yi, zi; 1/2, 22) for zi — > +0. 
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